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摘 要 : 柴 达 木 贫 地 内 的 湖泊 对 维系 区 域 生态 平衡 满足 生产 及 生活 用 水 ,保护 生态 环境 具有 重要 的 作用 。 随 着 
区 域 气候 变化 和 人 类 活动 影响 的 增强 ,盆地 内 湖泊 格局 发 生 了 一 系列 变化 ,为 了 分 析 湖 泊 变 化 特征 , 查 明 变化 原 
,以 快速 ,精确 .适用 性 强 为 目标 ,基于 Landsat 系 列 遥 感 影像 ,通过 分 析 可 鲁 克 湖 流域 湖泊 水 体 TOA(Top-of-atmo- 


sphere) 反 射 率 与 其 他 地 物 TOA 反射 率 的 差异 ,提出 了 一 种 湖泊 水 体 自动 提取 方法 一 一 湖泊 水 体 差分 模型 。 利 用 该 


模型 处 理 了 柴 达 木 盆地 的 近 百 个 Landsat 图像 场景 ,提取 了 不 同时 间 节 点 ,不同 空间 位 置 的 湖泊 水 体 ,使 用 总 体 分 类 
精度 Kappa 系数 及 用 户 精 度 对 其 精度 进行 评价 ,并 与 NDWI(Normalized Difference Water Index ) .MNDWI(Modified 
Normalized Difference Water Index) 方 法 对 该 地 区 湖泊 水 体 提取 结果 进行 了 比较 分 析 。 结 果 表明 :(1) 利用 TOA 反射 
率 差异 ,可 区 分 目标 与 非 目 标 地 物 ; (2) 基于 稳定 的 阔 值 ,湖泊 水 体 差 分 模型 可 实现 湖泊 水 体 信息 快速 提取 ,与 
NDWI、MNDWI 等 水 体 信息 提取 方法 相 比 ,能 够 更 加 有 效 抑制 地 表 河 流 、 冰 雪 、 阴 影 ,沼泽 湿地 等 干扰 因素 ,在 模型 
应 用 的 区 域内 ,平均 总 体 分 类 精度 与 用 户 精 度 均 达 到 99% 以 上 ,Kappa 达 到 0.9877;(3) 湖泊 水 体 差分 模型 的 输入 数 
据 既 可 以 是 TOA 反射 率 , 也 可 以 是 Landsat(Level-2) 的 地 表 反 射 率 数 据 ;(4) 湖泊 水 体 差 分 模型 适用 于 柴 达 木 盆地 
内 大 范围 区 域 的 湖泊 水 体 提 取 , 能 够 为 湖泊 水 体 动态 变化 规律 研究 提供 技术 支持 。 

关键 词 : Landsat; 湖泊 ; TOA 反射 率 ; 湖泊 水 体 差分 模型 ， 柴 达 木 盆地 


湖泊 参与 流域 水 循环 ,对 维持 区 域 水 量 平衡 ""、 
满足 生产 生活 用 水 、 河 流 流量 调节 ”发挥 着 重要 
作用 。 烷 达 木 分 地 地 处 青藏 高 原 东北 部 , 例 地 内 湖 
泊 众 多 , 且 多 位 于 盆地 中 心 的 低洼 地 带 ”, 为 地 表 与 
地 下 水 的 汇集 区 和 河流 水 系 的 尾 关 ,是 气候 变化 与 
人 类 活动 的 敏感 指示 顺 ” ,湖泊 及 其 引起 的 生态 环 
境 的 变化 在 一 定 程度 上 可 以 反映 区 域 乃 至 全 球 的 
气候 事件 ” ,在 干旱 的 内 陆 盆地 区 ,湖泊 是 维持 区 域 
生态 系统 的 重要 文 撑 ,同时 ,盆地 内 的 盐湖 作为 资 
源 也 是 当地 的 经 济 文 柱 ”"。 近 年 来 , 随 着 全 球 气候 
变化 和 人 类 活动 干扰 的 加 剧 , 侈 地 内 湖泊 出 现 了 面 
积 萎 缩 乃 至 干 润 的 现象 *”。 湖 泊 面 积 的 变化 不 仅 
影响 当地 的 社会 经 济 和 生态 环境 的 可 持续 发 展 ， 
而 且 对 区 域 乃 至 青藏 高 原 气 候 变 化 均 具有 重要 影 
Me] °° 。 因 此 ,精准 监测 上 某 达 木 盆地 湖泊 水 体 动态 ， 
揭示 自然 与 人 类 活动 对 湖泊 的 影响 规律 ,对 合理 开 
发 利用 和 保护 湖泊 水 域 ,维系 区 域 生态 平衡 具有 重 
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要 的 意义 "1。 

卫星 遥感 技术 以 其 探测 范围 广 、 信 息 容 量 大 、 
实时 性 与 动态 性 强 .免费 共享 等 优势 ,在 湖泊 水 体 
言 息 提取 中 得 到 了 国内 外 学 者 的 青睐 天 ,湖泊 水 
体 的 提取 最 早 主要 依靠 目 视 解 译 , 目 视 解 译 提 取水 
体 的 效果 好 坏 受 人 为 主观 因素 影响 ,而 且 耗 费 的 时 
间 成 本 ,人 力 成 本 高 。 用 计算 机 自动 提取 水 体 的 早 
期 研究 主要 是 利用 水 体 在 单个 波段 上 表现 出 的 与 
非 水 体 截 然 不 同 的 特性 ,如 1985 年 ,Jupp 等 "通过 
对 TM7(TM 影像 的 第 七 波段 ) 设 定 阔 值 识 别 水体 信 
息 ;1992 年 陆 家 驹 等 "通过 对 TM5 设 定 阅 值 识别 水 
体 信息 , 受 限 于 分 辨 率 , 对 于 1000~4000 m Z fal 1 
水 体 识别 起 来 还 有 困难 。 随 后 的 研究 利用 了 水 体 
在 多 个 波段 上 的 组 合 特征 ,产生 了 多 种 水 体 指数 方 
法 ,以 1996 年 McFeeters" 提 出 的 归 一 化 差异 水 体 
指数 NDWI(Normalized Difference Water Index) Aft 
表 , 通 过 抑制 植被 和 土壤 信息 ,提取 了 水 体 信 息 , 但 
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其 对 建筑 物 .阴影 r 细小 河流 的 区 分 度 不 高 ; 
杨 存 建 等 "利用 水 体 在 TM2 波 段 与 TM3 波 段 的 亮 
度 值 之 和 大 于 TM4 波段 与 TM5 波段 的 亮度 值 之 和 
这 一 特征 提取 水 体 ;人 徐 涵 秋 提出 了 改进 的 归 一 化 
差异 水 体 指数 法 MNDWI(Modified NDWI) , 较 好 地 
提取 了 城镇 水 体 信 息 ,该 方法 仅 适用 于 拥有 SWIR 
的 传 感 絮 , 易 受 冰雪 的 影响 ; 辣 需 等 ”提出 增强 型 
水 体 指数 EWI(Enhanced Water Index) ,在 提取 水 体 
时 区 分 了 半 干 润 河道 与 背景 噪声 ,但 对 细小 水 体 分 
离 效果 差 ,出现 湖泊 漏 提现 象 ; 丁 凤 ”将 TM7 波 段 
应 用 于 构建 新 型 水 体 指 数 NWI(New Water Index) 
上 ,能 够 清晰 提取 水 体 信息 ,但 无 法 区 分 山体 ;2014 
年 王 晴 晴 等 ?提出 简单 比值 型 水 体 指数 SRWI(Sim- 
ple Ratio of Water Index ) 用 来 增强 水 体 与 植被 .建筑 
物 、 裸 地 和 阴影 的 差异 ; 丁 占 峰 等 ”基于 对 避 阳 湖 
水 体 信息 提取 的 研究 提出 ONDWI( Override NDWI) 
用 来 去 除 水 体 中 含有 的 沼泽 湿地 ,同时 抑制 植被 和 
土壤 信息 ,但 对 冰雪 和 阴影 的 区 分 度 不 高 。 此 外 ， 
在 沙漠 水 体 信息 提取 ”城市 水 体 信息 提 取 *””、 
水 域 边界 或 河流 信息 提取 "背景 噪音 去 除 汪 
及 提高 土壤 .植被 与 水 体 信息 区 分 度 ”” 等 方面 ,学 
者 们 也 做 了 较为 深入 的 研究 ,建立 了 一 系列 有 效 实 
用 的 水 体 指 数 。 针 对 柴 达 木 贫 地 湖泊 水 体 动态 变 
化 ,前 人 主要 利用 目 视 解 译 . 单 波段 国 值 法 NDWTI、 
MNDWI 等 方法 开展 研究 ,取得 了 一 系列 成 果 ”“*。 
以 上 水 体 提 取 方 法 中 , 单 波段 国 值 法 与 比值 法 依赖 
动态 的 阔 值 对 湖泊 水 体 信息 进行 提取 , 且 对 阴影 、 
冰雪 噪声 的 去 除 效 果 不 理想 ;基于 各 种 水 体 指数 法 
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提取 水 体 信 息 的 羡 值 基本 稳定 为 正 值 ,因此 将 经 水 
体 模型 运算 后 结果 大 于 0 的 区 域 判 别 为 水 体 来 提取 
言 息 ,然而 试验 表明 NDWI 在 提取 西部 湖泊 水 体 信 
息 时 会 将 沼泽 湿地 与 冰雪 一 起 提取 出 来 ,MNDWI、 
EWI、NWI、ONDWI 等 方法 则 存在 湖泊 漏 提 和 误 提 
冰雪 的 情况 。 

随 着 信息 处 理 .定量 遥感 等 相关 技术 的 发 展 ， 
遥感 解 译 正 由 半自动 化 向 全 自动 智能 化 解 译 方向 
发 展 扩 ,对 遥感 解 译 的 速度 、 分 类 精度 与 适用 性 提 
出 了 更 高 的 要 求 。 为 了 解决 上 述 问题 ,本文 以 青海 
省 德 令 哈市 境内 的 可 鲁 克 湖 流域 为 试验 区 ,以 可 和 鲁 
克 湖 、 托 素 湖 、 厌 海 为 试验 对 象 建立 水 体 识别 模型 ， 
为 了 验证 模型 的 识别 精度 ,选择 了 柴 达 木 盆 地 中 的 
大 柴 晶 湖 .小 柴 旦 湖区 斯 库 勒 湖 和 察 尔 汗 盐湖 4 个 
典型 湖泊 作为 模型 的 验证 区 域 。 同 时 ,为 了 验证 模 
型 在 柴 达 木 盆 地 以 外 区 域 的 适用 性 ,将 该 模型 应 用 
于 邻近 柴 达 木 盆 地 的 青海 湖 ZTE, WA R 
池 和 法 池 的 水 体 提 取 中 。 基 于 Landsat 系列 卫星 影 
像 ,进行 湖泊 水 体 信 息 自动 提取 方法 研究 ,以 精度 
高 .速度 快 ` 适 用 性 强 为 目标 , 旨 在 提出 一 种 适合 柴 
达 木 盆地 区 域 的 自动 提取 湖泊 水 体 信 息 的 方法 , 进 
而 为 湖泊 水 体 的 动态 监测 提供 可 靠 的 技术 支撑 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 
柴 达 木 合 地 位 于 青藏 高 原 东 北 缘 , 地 处 青海 省 
西北 部 (34°40'~39°20'N .90°00'~99°20 下 ,图 1)。 中 
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图 1 研究 区 位 置 及 采样 点 分 布 
Fig. 1 Location of the study area and the distribution of sampling points 
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部 的 平原 区 海拔 2676~3000 m, 地 势 低洼 ,四 周 高 山 
环绕 ,西南 暖 湿 气流 难以 进入 ,降水 稀少 蒸发 量 
大 日照 充 足 , 属 于 内 陆 极端 干旱 高 寒气 候 。 盆 地 
内 共有 大 小 河流 79 条 ,这 些 河流 以 盆地 为 中 心 , 呈 
聚合 状 分 布 ,河水 通过 地 表 及 地 下 径流 最 终 注 入 各 
自 的 汇 水 中 心 ,形成 尾 闻 湖 ” 。 

1.2 数据 来 源 与 方法 

1.2.1 RAH 选用 Landsat 系 列 卫 星 影 像 数据 
作为 数据 源 ,数据 从 美国 地 质 勘探 局 网 站 (https:/ 
earthexplorer.usgs.gov/) 下载, 空间 分 辨 率 30 m ,影像 
级 别 为 Level 1TP( 地 形 校正 ) 级 。 除 了 柴 达 木 倪 地 
内 的 影像 外 ,还 收集 了 柴 达 木 盆地 内 外 其 他 区 域 不 
同时 期 ,不同 传 感 需 的 共计 90 个 Landsat 图 像 场 景 
数据 , 旨 在 验证 本 文 提 出 的 湖泊 水 体 提 取 方 法 在 不 
同时 间 ,空间 以 及 传感器 的 善 志 适用 性 。 本 文选 取 
自 2000 一 2020 年 可 鲁 克 湖 .小 上 某 旦 湖 ZREN 
察 尔 汗 盐湖 .青海 湖 等 区 域 的 部 分 影像 进行 论述 ， 
数据 基本 信息 见 表 1。 

1.2.2 还 感 数据 预 处 理 h TREH DN (digital 
number) 是 由 辐射 亮度 值 经 过 线性 变换 拉 伸 到 (0~ 
255) 范 围 内 后 得 到 的 灰 度 值 , 不 能 准确 表征 地 物 在 
每 个 波段 真实 的 反映 率 ,因此 在 ENVI 软 件 中 利用 
产品 提供 的 影像 定 标 文件 将 DN 值 转换 成 TOA(Top- 
of-atmosphere) 反 射 率 。TOA 反射 率 相对 于 DN 值 有 
两 大 优点 :一 是 消除 了 传感器 在 不 同时 间 获 取 数 据 
时 太阳 天 顶 角 带 来 的 余弦 效应 ;二 是 补偿 了 外 大 气 
太阳 辐 照 度 的 不 同 值 ”。 

1.2.3 采样 点 布设 为 了 获取 不 同 地 物 的 波谱 曲 
线 , 在 对 影像 进行 图 像 增 强 的 基础 上 ,结合 谷歌 地 
球 上 同 区 域 高 分 辩 率 影像 ,采用 目 视 判读 法 共 收 集 
了 55615 个 样本 点 ,其 中 包括 26676 个 水 体 样 点 ， 
28939 个 非 水 体 样 点 。 


此 外 ,在 青海 湖区 域 2018 年 7 月 的 影像 发 布 
后 ,于 当月 在 距 青 海 湖 东北 3 km 的 区 海 湖 周 边 , 对 
多 海 湖 水 体 和 非 水 体 之 间 的 边界 线 采 集 了 8406 个 
实测 定位 点 数据 ,这 些 采 样 点 组 成 的 虚线 代表 实际 
的 湖岸 线 ,用 于 湖泊 水 体 信息 提取 结果 的 视觉 评估 
以 及 精度 评价 参考 。 
1.2.4 技术 路 线 模型 建立 分 三 步 完成 :第 一 步 , 基 
于 下 载 的 Landsat 系列 影像 ,借助 ENVI 软 件 ,利用 
Landsat 产品 提供 的 元 数据 文件 ,将 DN 值 转换 为 
TOA 反射 率 值 ;第 二 步 , 通 过 实地 调研 ,结合 目 视 解 
译 ,构建 不 同 地 物 的 波谱 曲线 ,分 析 湖 泊 水 体 与 干 
扰 地 物 要 素 之 间 的 差异 ,初步 提出 湖泊 水 体 信息 提 
取 模 型 ;第 三 步 , 对 模型 提取 精度 进行 验证 ,检验 精 
度 超过 95% 时 ,确定 建立 湖泊 水 体 信息 提取 模型 ， 
检验 精度 较 低 时 ,重复 第 二 和 第 三 步 操 作 , 湖 泊 水 
体 信息 提取 模型 建立 的 思路 见 图 2。 
1.2.5 精度 评价 方法 ”为 了 验证 湖泊 水 体 信 息 提 取 
的 精度 ,引入 视觉 评估 法 和 混 清和 矩阵 ”2 种 方法 。 

混 痛 和 矩阵“” 中 的 总 体 分 类 精度 OA( 公 式 1)， 
Kappa 系数 (公式 2) ,以 及 用 户 精度 User accuracy 
(公式 3) 等 参数 可 以 定量 反映 分 类 效果 。 

oa = 正确 分 类 的 像 元 数 sop 


2a (1) 
参考 像 元 总 数 
q q 
ny ry = > Pata 
Kappa = 7 k=l (2) 


n ~ yi 
k=1 
RP: n EARMA; nj 是 第 类 地 物 样 本 点 
的 总 数 ; n, 是 从 分 类 数据 中 导出 第 类 样 点 的 总 
Bl; nm, 是 从 数据 集 导出 的 第 类 地 物 样本 点 的 总 
数 ; g 是 样本 类 别 的 总 数 。 


(3) 


n, 
User accuracy = — x 100 
Nis 


表 1 所 用 部 分 遥感 影像 基本 信息 


Tab.1 Basic information of some remote sensing images used 


地 区 feet 了 期 /年 -月 -日 条 带 号 行 号 云 量 /% 

可 重 克 湖 T™ 2000-04-10 135 34 13.45 
可 和 鲁 克 湖 ETM+ 2002-04-08 135 34 3 

可 和 鲁 克 湖 OLI 2016-07-27 135 34 2.36 

小 柴 旦 湖 TM 2004-06-15 136 34 0.79 
FR EAT) TM 2002-05-14 139 34 14.42 
青海 油 OLI 2018-07-19 133 34 20.49 
察 尔 汗 盐湖 OLI 2020-04-08 136 34 6.88 
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Landsat 系 列 影像 


ENVI 影像 定 标 文件 


获取 TOA 值 
ý 
构建 不 同 地 物 的 波谱 曲线 


曲线 是 否 存 在 差异 


构建 湖泊 水 体 信息 提取 模型 
检验 模型 精度 
y 


建立 湖泊 水 体 信息 提取 模型 


图 2 湖泊 水 体 信息 提取 模型 建立 的 技术 路 线 
Fig. 2 Technical route for establishing lake water body 


information extraction model 


总 体 分 类 精度 (OA) Kappa 系数 和 用 户 精度 的 
值 越 大 ,说 明 分 类 效果 越 好 。 


2 湖泊 水 体 提取 模型 的 提出 


2.1 湖泊 水 体 的 波谱 特征 

不 同类 型 的 地 物 对 一 定 波长 的 电磁 波 表现 出 
不 同 的 反射 特性 ,地 物 反射 电磁 波 特 性 在 遥感 影像 
上 的 定量 表现 就 是 TOA 反 射 率 数 值 不 同 ,通常 用 波 
谱 曲 线 来 刻画 地 物 在 不 同 波 段 的 TOA 反 射 率 表 
现 。 根 据 不 同 地 物 TOA 反射 率 的 差异 构建 模型 , 即 
通过 增强 目标 地 物 抑 制 非 目 标 地 物 的 方式 ,将 目标 
地 物 与 其 他 地 物 区 分 开 。 

在 上 某 达 木 贫 地 内 ,湖泊 水 体 的 矿 化 度 存 在 较 大 
的 差异 ,有 淡水 湖 . 半 咸 水 湖 咸水湖 和 盐湖 4 类 , 平 
均 高 程 约 2788 m”。 在 可 和 鲁 克 湖 流域 ,有 可 重 克 湖 
(淡水 湖 ) . 托 素 湖 (成 水 湖 ) 与 条 海 (盐湖 )3 个 湖泊 ， 
矿 化 度 分 别 为 :0.80~0.88 g:L'、29.62~30.20 gL, 
103.6~156 gL ', 流 域 平均 高 程 2980 m, JEM 
地 形 地 貌 还 是 矿 化 度 上 看 ,可 重 克 湖 流域 地 物 要 素 
的 光谱 特征 在 柴 达 木 盆地 均 具 有 典型 性 。 为 此 ,本 
文选 择 可 和 鲁 克 湖 流 域 作 为 试验 区 开展 湖泊 水 体 信 
息 提 取 方 法 研究 。 
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基于 上 述 思 路 ,在 可 重 克 湖 周 边区 域 , 在 实地 
调研 的 基础 上 ,通过 目 视 解 译 选 择 对 湖泊 水 体 信息 
提取 干扰 较 大 的 沼泽 湿地 冰雪、 阴影 . 裸 地 、 建 筑 
物 等 5 类 典型 要 素 ,通过 样本 点 (图 1) 获 取 其 在 不 同 
波段 上 的 TOA 反 射 率 , 分 别 求 取 其 平均 TOA 反射 
率 ,以 波段 为 横 坐 标 ,TOA 反射 率 为 纵 坐标 ,绘制 5 
类 典型 地 物 和 湖泊 水 体 的 波谱 曲线 (图 3)。 


一 阴影 一 湖泊 一 雪 
一 沼泽 湿地 一 裸 地 一 建筑 物 


TOA 反 射 率 ( 无 量 纲 ) 
二 
aN 


0 1 1 I + 
Blue Green Red NIR 


波段 


SWIR-1 SWIR-2 


注 :Blue 为 蓝光 ;Green 为 绿 光 ;Red 为 红 光 ;NIR 为 近 红 外 ;SWIR1 
为 短波 红外 1;SWIR2 为 短波 红外 2。 
图 3 研究 区 6 种 地 类 的 波谱 曲线 


Fig.3 Spectral signatures of six major land use 


classes in the study area 


从 图 3 可 以 看 出 ,湖泊 水 体 的 TOA 反射 率 与 其 
他 5 类 地 物 有 明显 的 不 同 ,表现 在 以 下 2 个 方面 :一 
是 湖泊 是 唯 类 波谱 曲线 整体 呈 下 降 趋势 的 地 
物 ; 二 是 湖泊 在 Blue .Green 波段 的 反射 率 值 接近 其 
他 地 物 ,而 在 Red NIR、SWIR1 .SWIR2 波段 明显 低 
于 其 他 地 物 。 基 于 上 述 差异 ,在 研究 区 内 ,可 以 根 
据 TOA 反 射 率 来 区 分 湖泊 水 体 和 其 他 地 物 。 
2.2 湖泊 水 体 差 分 模型 的 提出 

对 于 清水 ,在 Blue-Green 波 段 反 射 率 为 4%~5%， 
波长 0.6 pm 以 上 的 Red 波 段 反 射 率 降 到 2%~3% ,在 
NIR 波 段 上 几乎 吸收 全 部 的 入 射 能 量 和 外 。 对 于 柴 达 
木 盆 地 的 湖水 来 说 , 随 着 悬浮 物 浓度 的 增加 ,与 清 
水 相 比 ,反射 率 数 值 会 发 生 一 定 的 影响 ” ,但 整体 
的 反射 率 峰 值 仍 然 处 于 Blue .Green 波段 范围 内 , 且 
反射 率 随 波 长 的 增加 而 减 小 “”。 基 于 此 ,对 1.2.3 
小 节 中 的 采样 结果 , 取 TOA 反射 率 ,执行 (Blue+ 
Green-Red-NIR-SWIR1-SWIR2) 运 算 ,结果 表明 ， 
湖泊 水 体 中 ,所 有 湖泊 样本 点 的 TOA 反射 率 运 算 结 
果 均 大 于 0, 而 几乎 所 有 的 其 他 地 物 样 本 点 表现 出 
相反 的 运算 结果 。 由 此 可 以 初步 得 出 以 下 结论 : 
(Blue+Green-Red-NIR-SWIR1-SWIR2) 运 算 可 以 
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将 湖泊 水 体 信息 增强 为 大 于 0 的 值 ,将 其 他 非 水 体 
信息 抑制 为 小 于 0 的 值 。 据 此 ,提出 了 基于 Landsat 
系列 卫星 影像 的 湖泊 水 体 信 息 提取 模型 一 一 湖泊 
水 体 差 分 模型 (Lake Water Differential Model, 
LWDM). 
Blue + Green — Red - NIR - SWIR1 - SWIR2 > 0 
(4) 


3 ”结果 与 分 析 


为 了 验证 模型 的 可 行 性 ,对 比 不 同 模型 的 提取 
效果 ,本 文选 择 位 于 柴 达 木 贫 地 东部 的 可 鲁 克 湖 、 
条 海 和 托 素 湖 ,中 部 的 大 上 某 旦 湖 .小柴 旦 湖 ,西部 的 


(a) 真 彩色 合成 影像 
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(c) NDWI 方 法 水 体 提 取 结 果 
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同 传感器 数据 ,对 LWDM 模型 进行 验证 。 同 时 ,为 
证 明 模型 在 更 大 空间 尺度 的 普遍 适用 性 ,还 选择 了 
临近 柴 达 木 伟 地 东部 边缘 的 青海 湖 、 帮 海 湖 作为 验 
证 对 象 。 
3.1 湖泊 水 体 信息 提取 结果 

采用 LWDM NDWI 和 MNDWT 方 法 部 分 提取 结 
果 如 下 图 所 示 ( 图 4~ 图 5)。 图 4a 采 用 真 彩色 影像 
合成 并 进行 平方 根 拉 伸 , 图 5a 采 用 真 彩色 影像 合成 
并 进行 高 斯 增强 ,由 在 突出 湖泊 在 原 影像 上 的 
位 置 。 

在 可 和 鲁 克 湖 流 域内 ,由 图 4 可 知 ,3 种 方法 均 可 


(b) LWDM 方 法 水 体 提取 结果 
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Fig. 4 Results of lake water extraction in Keluke Lake 
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(a) 真 彩色 合成 影像 
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图 5 青海 湖区 域 湖泊 水 体 提取 结果 
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(b)LWDM 方 法 水 体 提取 结果 
100°30'0"E 


100°35'0"E 
Bee fas, 


ZRH 


37°0'0"N 


”野外 实测 海岸 线 
ne 


= 
100°30'0"E 100°35'0"E 
(d) MNDWI 方 法 水 体 提 取 结 果 
100°30'0"E 100°35'0"E 
T a 


ARE 


37°0'0"N 


We 
一 野外 实测 海岸 线 
一 提取 结果 


100°30'0"E 100°35'0"E 


Fig. 5 Results of lake water extraction in Qinghai Lake 


以 将 湖泊 水 体 信 息 提 取出 来 ,但 NDWI 方 法 在 提取 
湖泊 信息 (图 4c) 的 同时 提取 了 河流 水 体 、 沼 泽 湿地 
和 部 分 阴影 的 信息 ,MNDWI 方 法 同时 提取 了 河流 
水 体 、 沼 泽 湿 地 、 阴 影 和 冰雪 的 信息 (图 4q); 这 种 情 
况 在 小 柴 且 湖 周边 区 域 以 及 系 斯 库 勒 湖区 域 同 样 
存在 。 

在 青海 湖区 域 , 选 择 同 一 时 间 段 (2018 年 7 月 )， 
对 比 野 外 实测 湖岸 线 与 基于 LWDM 的 湖泊 提取 结 
果 ( 图 6), 可 以 看 出 ,提取 结果 与 实测 湖岸 线 误差 在 
一 个 像素 内 ,能 够 反映 实际 湖水 岸 线 的 曲折 情况 。 

里 然 , 在 模型 建立 之 初 就 已 经 考虑 到 了 烷 达 木 
盆地 不 同 湖泊 水 体 矿 化 度 的 差异 ,作为 测试 对 象 的 
录 海 ,其 矿 化 度 已 经 达到 了 元 水 的 标准 ,但 与 察 尔 
汗 盐 湖 矿 化 度 高 达 400 g*17' 相 比 ,还 存在 较 大 的 差 
BE ,为 了 进一步 验证 模型 的 适用 性 ,运用 本 文 提出 
的 方法 ,提取 了 察 尔 汗 盐湖 及 其 周边 的 湖泊 水 体 信 
息 ( 图 7) ,由 网 7 可 知 ,LWDM 较为 完整 的 提取 了 湖 
泊 水 体 信息 ,有 效 去 除了 区 内 的 道路 、 建 筑 物 等 干 


扰 信 息 , 总体 分 类 精度 为 97.98% , Kappa 系数 为 
0.9339, 用 户 精度 99.94% ,说 明 LWDM 模型 在 柴 达 
木 盆地 内 具有 和 较 好 的 适用 性 。 

通过 上 述 对 比 可 以 得 出 如 下 结论 :LWDM 在 提 
取 湖 泊 水 体 信息 的 同时 ,不 仅 可 以 去 除 植被 . 裸 地 、 
建筑 物 等 干扰 信息 ,也 可 以 较 好 地 去 除 绝 大 部 分 
雪 阴影、 沼泽 湿地 和 地 表 河 流 等 干扰 因素 ,与 
NDWI, MNDWI 等 水 体 提 取 方 法 相 比 ,LWDM 方法 
可 以 更 好 地 提取 柴 达 木 合 地 内 的 湖泊 水 体 信 息 。 

3.2 精度 评价 

根据 原始 的 混淆 矩阵 ,经 过 公式 (1) ~ (3) 的 运 
算 ,得 到 精度 评价 结果 ( 表 2)。 

由 表 2 可 知 ,湖泊 水 体 信息 的 提取 精度 表现 为 : 
LWDM>NDWI>MNDWI。 在 青海 湖区 域 附 近 ,3 种 
方法 提取 湖泊 水 体 信 息 的 精度 接近 , 而 在 柴 达 木 盆 
地 内 的 可 重 克 湖区 域 .小 柴 旦 湖区 域 . 胡 斯 库 勒 湖 
区 域 . 察 尔 汗 盐湖 区 域 ,3 种 水 体 提取 方法 在 湖泊 提 
取 精 度 上 表现 出 显著 的 差异 ,其 原因 是 NDWT、 
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图 6 青海 湖区 域 LWDM 方 法 水 体 局 部 提取 结果 
Fig.6 Results of lake water extraction by LWDM in Qinghai Lake 
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(b) LWDM 方 法 水 体 提取 结果 
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图 7 察 尔 汗 盐湖 水 体 提取 结果 
Fig.7 Results of lake water extraction by LWDM in Chaerhan Salt Lake 


MNDWI 在 提取 湖泊 水 体 时 ,会 同时 提取 地 表 河 流 、 
冰雪 沼泽 湿地 等 非 湖泊 水 体 信 息 , 从 而 降低 了 用 
户 精度 。 

此 外 ,在 柴 达 木 盆地 的 其 他 区 域 不 同时 期 .不 
同 传感器 共计 90 个 Landsat 图 像 场景 数据 中 ,OA 达 
到 99.48% ,Kappa 值 达到 0.9877。 

由 此 可 以 得 出 如 下 结论 :基于 Landsat 系 列 影 
像 ,LWDM 可 以 较 好 地 提取 柴 达 木 岔 地 湖泊 水 体 信 
息 , 相 对 于 NDWI.MNDWI 等 水 体 信息 提取 方法 ,去 


噪 能 力 强 ,提取 精度 显著 提高 。 对 于 上 某 达 木 盆地 之 
外 的 其 他 区 域 , 针 对 不 同 清洁 度 的 湖泊 ,LWDM 方 
法 的 适用 性 有 待 进一步 的 研究 。 
3.3 提取 结果 分 析 

以 Landsat TM ETM+ 和 0OLI 多 光谱 遥感 影像 为 
数据 源 ,通过 分 析 可 和 鲁 克 湖 流域 内 湖泊 水 体 与 冰 
雪 阴影 沼泽 湿地 、 裸 地 、 建 筑 物 等 对 象 的 TOA 反 
射 率 差 异 ,根据 湖泊 水 体 在 Blue .Green 波段 的 TOA 
反射 率 高 于 其 余波 段 , 且 不 同 波段 TOA 反射 率 值 呈 
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表 2 湖泊 水 体 信息 提取 精度 评价 结果 


Tab.2 Evaluation result of information extraction accuracy of lake water 


试验 区 域 传感器 成 像 时 间 / 年 -月 -日 方法 OA/% Kappa 用 户 精度 /% 

可 重 克 湖 TM 2000-04-10 LWDM 99.65 0.9958 99.39 
可 鲁 克 湖 ETM+ 2002-04-08 99.43 0.9932 100 

可 重 克 湖 OLI 2016-07-27 99.87 0.9987 99.47 
小 业 旦 湖 TM 2004-06-15 99.65 0.9948 100 

ARS i TM 2002-05-14 99.91 0.9988 99.65 
青海 油 OLI 2018-07-19 99.87 0.9984 99.20 
察 尔 汗 盐湖 OLI 2020-04-08 97.98 0.9339 99.94 
平均 值 99.48 0.9877 99.66 
可 重 克 党 TM 2000-04-10 NDWI 95.42 0.9450 79.12 
DE re ETM+ 2002-04-08 95.07 0.9409 77.18 
可 重 克 党 OLI 2016-07-27 97.83 0.9728 90.20 
58 Ait TM 2004-06-15 95.60 0.9340 88.34 
IRIT PEA TM 2002-05-14 96.38 0.9517 87.34 
青海 油 OLI 2018-07-19 99.63 0.9956 97.84 
察 尔 汗 盐湖 OLI 2020-04-08 96.27 0.7363 99.89 
平均 值 96.60 0.9252 88.56 
可 重 克 党 TM 2000-04-10 MNDWI 77.67 0.7320 42.74 
ETEW ETM+ 2002-04-08 64.77 0.5773 32.12 
可 重 克 党 OLI 2016-07-27 96.80 0.9599 86.19 
小 柴 量 济 TM 2004-06-15 65.74 0.4861 49.31 
IRIT ETD TM 2002-05-14 71.10 0.6147 46.38 
青海 油 OLI 2018-07-19 91.10 0.8932 65.19 
察 尔 汗 盐湖 OLI 2020-04-08 96.27 0.7301 89.91 
平均 值 80.49 0.7133 58.83 


整体 下 降 趋 势 的 特征 ,利用 蓝光 ARG LG ZT. 
外 与 短波 红外 波段 ,构建 了 柴 达 木 盆地 湖泊 水 体 信 
息 提取 模型 一 一 湖泊 水 体 差分 模型 (LWDM)。 将 该 
模型 应 用 于 柴 达 木 倪 地 的 可 和 鲁 克 湖 、 托 素 湖 、 胡 海 、 
阴 斯 库 勒 湖 、 小 柴 旦 湖 以 及 盆地 附近 青海 湖 地 区 的 
大 海 湖 等 湖泊 水 体 信 息 的 提取 ,与 Mcfeeters 的 
NDWI 和 人 徐 涵 秋 ”的 MNDWI 相 比 ,LWDM 不 仅 可 以 
较为 准确 地 提取 出 湖泊 水 体 信息 ,而 且 可 以 有 效 去 
除 冰雪 、 了 明 影 ,沼泽 湿地 、 建 筑 物 等 干扰 信息 。 模 型 
提出 后 ,利用 Landsat(Level-2 产品 ) 的 地 表 反 射 率 
作为 输入 数据 对 模型 进行 进一步 测试 ,结果 表明 其 
提取 精度 与 TOA 反射 率 作 为 输入 数据 的 提取 精度 
基本 一 致 。 因 此 ,本 文 提出 的 湖泊 水 体 差 分 模型 既 
可 以 将 TOA 反 射 率 作为 输入 数据 ,也 可 以 将 地 表 反 
射 率 作为 输入 数据 ;同时 应 用 不 同时 段 不同 季节 
的 Landsat 系列 影像 数据 进行 了 测试 ,结果 表明 ,该 
模型 所 反映 的 湖泊 水 体 规律 在 不 同年 份 ,不同 季 节 


注 : 加 粗 的 数据 表示 不 同 水 体 提取 方法 对 不 同 湖泊 水 体 信息 提取 结果 精度 的 平均 值 。 


是 一 致 的 。 此 外 ,NDWI 和 MNDWI 提 出 的 目的 是 提 
取 地 表 水 体 信息 ”” ,并 未 考虑 湖泊 水 体 信息 与 河 
流水 体 信 息 的 差异 ,因此 ,在 提取 湖泊 水 体 信息 的 
同时 ,也 混入 了 地 表 河 流水 体 信 息 ( 图 4c~ 图 4d)。 
而 LWDM 建立 的 目标 就 是 提取 湖泊 水 体 信 息 , 在 前 
期 试验 研究 过 程 中 充分 考虑 了 湖泊 水 体 与 地 表 水 
体 信 息 之 间 的 差异 ,在 提取 湖泊 水 体 信息 的 同时 ， 
较 好 地 剔除 地 表 河 流 的 水 体 信 息 ( 图 4b)。 

从 提取 精度 评价 的 结果 来 看 ,LWDM 方 法 的 提 
取 精 度 要 高 于 NDWI 和 MNDWI, 原 因 有 三 :第 一 ， 
NDWI 是 在 分 析 植 被 与 水 体 在 可 见 光 波段 和 近 红 外 
波段 的 反射 强度 特征 后 提出 的 ,目的 是 最 大 程度 
地 抑制 植被 信息 ,增强 水 体 信 息 。 而 MNDWI 是 针 
对 NDWI 不 能 很 好 区 分 城市 范围 内 的 水 体 而 改进 的 
一 种 水 体 指数 ” ,目的 是 增加 建筑 物 和 水 体 的 差 
异 , 进 而 更 精确 的 提取 城市 范围 内 的 水 体 ,但 是 不 
能 很 好 地 抑制 含水 量 较 大 的 土壤 干扰 信息 。 二 者 
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以 提取 地 表 水 体 (包括 河流 、 湖 泊 水 库 等 ) 为 目标 ， 
建立 时 更 多 地 关注 植被 ,建筑 物 对 水 体 信息 提取 的 
干扰 。 第 二 ,虽然 徐 涵 秋 在 提出 MNDWI 时 ,选择 
河流 .湖泊 海洋 3 种 水 体 对 NDWI 和 MNDWI 进 行 
了 水 体 信 息 提 取 试 验 ,实验 结果 证 明 ,MNDWI 具 有 
一 定 的 阴影 去 除 效 果 , 但 未 考虑 像 某 达 木 盆地 这 样 
复杂 的 环境 :一 是 使 地 属于 高 原型 内 陆 盆地 , 极 高 
山高 山 戈壁 .丘陵 平原. 盐 沼 等 地 貌 并 存 ;二 是 
区 内 湖泊 河流 密布 ,沼泽 湿地 普遍 存在 ,不 同 湖泊 
矿 化 度 差 异 较 大 ;三 是 盆地 周边 的 极 高 山 .高 山区 
有 冰雪 履 善 区 。 这 在 一 定 程度 上 降低 了 MNDWI 模 


利用 LWDM 对 其 周边 区 域 的 湖泊 水 体 信息 进行 了 
提取 ,借助 总 体 分 类 精度 OA Kappa 系数 及 用 户 精 
度 对 提取 结果 进行 了 精度 评价 ( 表 3), 由 表 3 可 知 ， 
天 池 的 提取 结果 较 好 ,太湖 和 滇池 的 效果 稍 差 , 原 
是 太湖 和 滇池 的 水 质 相对 较 差 ,说 明 LWDM 模 型 
对 清洁 水 体 的 适用 性 较 好 ,而 对 受 污染 水 体 的 适用 
性 相对 较 差 ,前 期 笔者 对 太湖 含 蓝藻 水 体 进 行 了 专 
题 研 究 ,提出 了 基于 Landsat 影 像 的 含 蓝 落水 体 提 取 
方法 (Blue+Green-Red-SWIR1-SWIR2>0)50 ,通过 
对 比 发 现 , 提 取 结 果 优 于 NDWI、MNDWI、MBWI、 
NWIWLos, 后 续 针 对 不 同 清洁 度 的 湖泊 水 体 信息 ， 


型 提取 柴 达 木 贫 地 湖泊 水 体 信 息 的 精度 。 第 三 ， 
LWDM 模 型 以 提取 柴 达 木 盆地 湖泊 水 体 信 息 为 目 
标 ,建立 时 不 仅 考虑 了 上 述 背 景 条 件 ,而 且 借助 波 
谱 曲 线 捕获 到 了 湖泊 水 体 与 其 他 地 物 要 素 的 TOA 
反射 率 差 异 , 模 型 在 精确 刻画 地 物 差异 的 同时 , 提 
高 了 提取 结果 的 可 靠 性 。 

LWDM 在 建立 时 既 考 虑 了 可 和 鲁 克 湖 流域 内 水 
体 在 可 见 光 波段 ,特别 是 蓝光 与 绿 光波 段 反射 率 相 
对 较 强 ,在 近 红 外 与 短波 红外 强 吸 收 的 特征 ,又 考 
虑 了 DN 值 在 表征 地 物 时 失真 的 不 足 。 模 型 建立 过 
程 中 采用 差 值 运算 方法 并 引入 较 多 的 波段 ,一 方面 
可 以 有 效 避 人 免 由 于 比值 运算 造成 局 部 信息 丢失 的 
现象 ; 另 一 方面 引入 多 波段 可 以 探测 到 更 多 的 水 体 
微细 变化 信息 ,对 与 柴 达 木 盆地 相似 的 不 同 矿 化 度 
湖泊 水 体 信息 并 存 的 区 域 ,可 取得 比 NDWI、MND- 
WI 更 好 的 提取 效果 。 在 净值 选择 方面 ,有 学 者 采用 
神经 网 络 技术 自动 计算 样本 阔 值 ,得 到 了 水 体 信 息 
自动 提取 模型 ,但 该 方法 容易 得 到 局 部 极 小 值 而 延 
长 计算 时 间 ,而 本 文 提 出 的 LWDM 模型 将 提取 湖泊 
水 体 信 息 的 浆 值 设 定 为 0, 与 NDWI 和 MNDWI 在 局 
部 地 区 信息 提取 时 需要 调整 闪 值 ”操作 及 采用 神 
经 网 络 技术 导致 计算 时 间 较 长 相 比 ,LWDM 提取 湖 
泊 水 体 信 息 操 作 更 加 便捷 。 

为 了 考证 LWDM 的 适用 性 ,选择 柴 达 木 盆 地 以 
外 位 于 不 同 地 区 的 天 池 太湖 滇池 作 为 验证 对 象 ， 


还 需要 对 模型 进行 完善 ,以 提高 其 适用 性 。 利 用 
MODIS 数据 (成 像 时 间 ,2005 年 5 月 25 日 ) 进 一 步 对 
LWDM 的 适用 性 进行 分 析 ,选择 盆地 内 外 的 哈 拉 
湖 .青海 湖 作为 试验 对 象 ,对 提出 的 水 体 提取 指数 
进行 了 试验 验证 ,采用 LWDM 提取 了 2 个 湖泊 水 体 ， 
提取 结果 的 总 体 分 类 精度 为 99.81% ,Kaapa 系数 为 
0.9973 , 用户 精度 99.78% ,说 明 对 于 不 同 的 传 感 需 ， 
企 青藏 高 原 地 区 该 方法 也 有 较 好 的 适用 性 。 

由 于 湖泊 水 体 信息 受制 于 多 种 因素 ,当前 的 研 
究 方 法 尚 不 能 完成 全 面 、 系 统 的 湖泊 水 体 信 息 提 
取 ,特别 是 对 较 差 水 质 的 湖泊 水 体 信 息 ,模型 还 有 
待 进一步 完善 。 本 文 基于 前 人 的 研究 成 果 , 对 柴 达 
木 盆地 的 湖泊 水 体 分 布 自 动 提取 方法 进行 了 探索 
性 研究 ,后 续 还 要 在 以 下 2 个 方面 加 强 研究 :一 是 引 
入 机 器 学 习 等 人 工 智能 技术 ,自动 识别 青藏 高 原 湖 
泊 水 体 与 冰雪 沼泽 湿地 、 阴 影 等 要 素 的 TOA 反 射 
率 差异 ,完善 湖泊 水 体 信 息 提取 模型 ,将 其 应 用 到 
塔里木 盆地 .云贵 高 原 .四 川 盆 地 等 不 同 区 域 , 在 应 
用 中 不 断 改进 模型 ,增强 模型 的 普 适 性 ;二 是 考虑 
SPOT , QuickBird , ASTER .IKONOS 及 国产 高 分 辩 率 
对 地 观测 等 卫星 遥感 数据 ,采用 多 源 遥 感 影像 数据 
融合 技术 ,关注 气象 变化 与 人 类 活动 的 影响 ,提高 
湖泊 水 体 分 布 信息 的 提取 精度 ,进而 为 湖泊 水 体 动 
态 变 化 及 其 周边 区 域 生 态 环 境 演 化 提供 理论 与 技 
术 支 撑 。 


表 3 滇池 ,太湖 .天池 湖 泊 水 体 信息 提取 精度 评价 结 


Tab. 3 Evaluation results of water body information extraction accuracy of Dianchi Lake, Taihu Lake and Tianchi Lake 


实验 区 域 传感器 成 像 时 间 / 年 -月 -日 
滇池 OLI 2021-03-25 
太湖 OLI 2021-06-23 
天 池 OLI 2021-06-18 


方法 OA/% Kappa JPE% 
LWDM 90.13 0.6295 89.86 
96.03 0.8185 96.32 
99.79 0.9413 99.82 
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AS SBE F Landsat 系列 遥感 影像 ,通过 对 比分 析 
可 和 鲁 克 湖 流 域 不 同 地 物 要 素 的 TOA 反射 率 差异 , 建 
立 了 基于 Landsat 影 像 的 湖泊 水 体 分 布 自动 提取 方 
法 一 一 湖泊 水 体 差分 模型 (LWDM) ,将 该 方法 应 用 

于 柴 达 木 盆地 TZ I te eE 
湖泊 水 体 信息 提取 ,并 分 析 了 同一 卫星 不 同 传 感 
器 不同 成 像 时 间 不 同 卫星 不 同 传感器 的 提取 结 
果 。 主 要 结论 如 下 : 

(1) 在 统计 分 析 的 基础 上 ,根据 地 物 在 遥感 影 
像 不 同 波段 的 TOA 反 射 率 差异 特征 ,通过 增强 目标 
地 物 信息 ,抑制 非 目标 地 物 信 息 , 可 以 有 效 提取 所 
需 信息 。 

(2) 基于 TOA 反射 率 差异 特征 提出 的 湖泊 水 体 

差分 模型 LWDM) 可 以 用 稳定 的 阔 值 (0) 实 现 区域 湖 
泊 水 体 信 aus 决 速 提 取 , 与 NDWI、MNDWI 等 方法 


相 比 ,在 柴 达 木 倪 地 内 ,该 方法 能 有 效 抑 制 地 表 河 
流 、 冰 雪 、 沼 泽 湿 地 、 胃 影 等 干扰 因素 ,显著 提高 湖 
泊 水 体 信 息 提取 的 精度 。 


(3) 基于 Landsat 系列 卫星 影像 ,LWDM 可 用 于 
柴 达 木 盆 地 湖泊 水 体 分 布 快速 识别 ,基于 不 同 传 感 
器 不 同时 期 的 遥感 影像 ,通过 对 比 不 同 湖泊 提取 结 
果 ,0A 达 到 99.48% ,Kappa 值 达到 0.9877 ,用 户 精度 
达到 99.66%。 基 于 不 同 卫星 传感器 ,利用 MODIS 32% 
感 影像 数据 ,提取 了 青海 湖 和 哈 拉 湖 水 体 信息 ,OA 
达到 99.81% , Kappa 值 达到 0.9973 ,用 户 精 度 达 到 
99.78% ,结果 具有 较 高 的 可 信和 度 ,可 为 柴 达 木 盆 地 
湖泊 水 体 动 态 监 测 及 湖泊 周边 生态 环境 保护 提供 
技术 文 持 。 
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An automatic method for delineating lake surfaces in 
Qaidam Basin using Landsat images 


WEN Guangchao, LI Xing, WU Bingjie, WANG Xiaohe, XIE Hongbo 


(Institute of Resource and Environment, Henan Polytechnic University, Jiaozuo 454000, Henan, China) 


Abstract: Lakes in the Qaidam Basin maintain regional ecological balance, supply production and domestic 
water consumption, and protect the ecological environment. These lakes have undergone a series of changes in 
response to enhanced regional climate change and human activities. To analyze characteristics of lake changes 
and investigate casual factors, an automatic method for lake surface delineation based on Landsat remote sensing 
images with a different series of sensors (the Lake Water Differential Model) was proposed by understanding 
TOA (top-of-atmosphere) reflectance differences between water and non-water surfaces in the Koruk Lake Basin. 
This model is applied to extract lake information at varying times or locations. Model performance is evaluated 
using overall accuracy, the Kappa coefficient, and user’s accuracy, and the derived lake surfaces were compared 
to those from the NDWI (Normalized Difference Water Index) and MNDWI (Modified Normalized Difference 
Water Index) methods. These results show the following: (1) target and non-target objects are distinguished by a 
TOA reflectance difference. (2) Based on a stable threshold, the Lake Water Differential Model can delineate lake 
surfaces. Compared with the NDWI, MNDWI, and other water body information extraction methods, the Lake 
Water Differential Model can more effectively suppress interference factors, such as surface rivers, ice, snow, 
shadows, swamps, and wetlands. In areas of model application, our proposed method can achieve a performance 
of 99.48% of the average OA, 99.66% of user accuracy, and 0.9877 of Kappa. (3) Input data of the model can be 
either TOA or Landsat (level- 2) surface reflectance data. (4) The model is suitable for extracting lake water 
information in a large area of Qaidam Basin and provides technical support for the study of lake surface dynamics. 


Keywords: Landsat; lake; TOA reflectance; lake water differential model; Qaidam Basin 


